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Abstract. We study a vacation-type queueing model, and a single-server multi-queue polling 
model, with the special feature of retrials. Just before the server arrives at a station there is 
some deterministic glue period. Customers (both new arrivals and retrials) arriving at the 
station during this glue period will be served during the visit of the server. Customers arriving 
in any other period leave Immediately and will retry after an exponentially distributed time. 
Our main focus is on queue length analysis, both at embedded time points (beginnings of 
glue periods, visit periods and switch- or vacation periods) and at arbitrary time points. 
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1 Introduction 


Queueing systems with retrials are characterized by the fact that arriving customers, who find the 
server busy, do not wait in an ordinary queue. Instead of that they go into an orbit, retrying to 
obtain service after a random amount of time. These systems have received considerable attention 
in the literature, see e.g. the book by Falin and Templeton na, and the more recent book by 
Artalejo and Gomez-Corral |3]. 

Polling systems are queueing models in which a single server, alternatingly, visits a number of 
queues in some prescribed order. Polling systems, too, have been extensively studied in the literature. 
For example, various different service disciplines (rules which describe the server’s behaviour while 
visiting a queue) and both models with and without switchover times have been considered. We 
refer to Takagi and Vishnevskii and Semenova |28| for some literature reviews and to Boon, 

van der Mei and Winands [B], Levy and Sidi m and Takagi |23| for overviews of the applicability 
of polling systems. 

In this paper, motivated by questions regarding the performance modelling of optical networks, 
we consider vacation and polling systems with retrials. Despite the enormous amount of literature on 
both types of models, there are hardly any papers having both the features of retrials of customers 
and of a single server polling a number of queues. In fact, the authors are only aware of a sequence 
of papers by Langaris |12I13I14| on this topic. In all these papers the author determines the mean 
number of retrial customers in the different stations. In m the author studies a model in which 
the server, upon polling a station, stays there for an exponential period of time and if a customer 
asks for service before this time expires, the customer is served and a new exponential stay period 
at the station begins. In m the author studies a model with two types of customers: primary 
customers and secondary customers. Primary customers are all customers present in the station at 
the instant the server polls the station. Secondary customers are customers who arrive during the 
sojourn time of the server in the station. The server, upon polling a station, first serves all the 
primary customers present and after that stays an exponential period of time to wait for and serve 
secondary customers. Finally, in mi the author considers a model with Markovian routing and 
stations that could be either of the type of m or of the type of m- 


This is an invited, considerably extended version of |10| . The main additions are Subsection 


2.4 


and 


Section|^ These present respectively the optimal behaviour of a single queue system, and the performance 
analysis for a general number of queues. 








In this paper we consider a polling station with retrials and so-called glue periods. Just before 
the server arrives at a station there is some deterministic glue period. Customers (both new arrivals 
and retrials) arriving at the station during this glue period "stick" and will be served during the 
visit of the server. Customers arriving in any other period leave immediately and will retry after 
an exponentially distributed time. 

Our study of queueing systems with retrials and glue periods was at first instance motivated by 
questions regarding the performance modelling and analysis of optical networks. Optical fibre offers 
some big advantages for communication w.r.t. copper cables: huge bandwidth, ultra-low losses, 
and an extra dimension - the wavelength of light. Performance analysis of optical networks is a 
challenging topic (see e.g. Maier [1^ and Rogiest (22]). In a telecommunication network, packets 
must be routed from source to destination, passing through a series of routers and switches. In 
copper-based transmission links, packets from different sources are time-multiplexed. This is often 
modeled by a single server polling system. In optical switches, too, one has the need for a protocol 
to decide which packet may be transmitted. One might again use a cyclic polling strategy, cyclic 
meaning that there is a fixed pattern for giving service to particular ports/stations. However, unlike 
electronics, buffering of optical packets is not easy, as photons cannot be stopped. Whenever there 
is a need to buffer photons, they are made to move locally in fiber loops. These fiber loops or fiber 
delay lines (FDL) originate and end at the head of a switch. When a photon arrives at the switch 
at a time it cannot be served, it is sent into an FDL, thereby incurring a small delay to its time of 
arrival without getting lost or displaced. Depending on the availability, requirement, traffic, size of 
photon and other such factors, the length (delay produced) of these FDLs can differ. Hence we 
assume that these FDLs delay the photons by a random amount of time. Also, if a packet does not 
receive service after a cycle through an FDL, then depending on the model it can go into either 
the same or a longer or a shorter or randomly to any of the available FDLs. Hence we assume that 
two consecutive retrials are independent of each other. This FDL feature can be modelled by a 
retrial queue. 

A sophisticated technology that one might try to add to this is varying the speed of light by 
changing the refractive index of the fiber loop, cf. m- By increasing the refractive index in a 
small part of the loop we can achieve ‘slow light’, which implies slowing the packets. When a 
port ‘knows’ that it will soon be served, it may start the process of increasing the refractive index 
at FDLs and at the end of fibers of incoming packets. By doing this, it slows down the packets 
which arrive at the station just before the visit period of the station begins. This feature is, in 
our model, incorporated as glue periods immediately before the visit period of the corresponding 
station. Packets arriving in this glue period can be served in that subsequent visit period. 

The concept of glue period is, to the best of our knowledge, new in polling systems. It may also 
be interpreted as a reservation period. We view a reservation period as a period in which customers 
can make a reservation at a station for service in the subsequent visit period of that station. In our 
case, such a reservation period occurs immediately before the visit period, and could be viewed as 
the last part, Gi, of a switchover period of length Si + Gi. Ordinary gated polling could be viewed 
as a service discipline in which it is always possible to make a reservation for the following visit 
period. 

The main contributions of the paper are the following, (i) For the case of a single queue with 
vacations and glue periods, we provide a detailed queue-length analysis at particular embedded 
epochs and at an arbitrary epoch. We also show how to choose the length of the glue period that 
minimizes the mean number of customers in the system, (ii) We also provide a detailed queue-length 
analysis for the Wqueue polling case - again at particular embedded epochs and at an arbitrary 
epoch. 

The paper is organized as follows. In Section we consider the case of a single queue with 
vacations and retrials; arrivals and retrials only "stick" during a glue period. We study this case 
separately because (i) it is of interest in its own right, (ii) it allows us to explain the analytic 
approach as well as the probabilistic meaning of the main components in considerable detail, (iii) it 
makes the analysis of the multi-queue case more accessible, and (iv) results for the one-queue case 
may serve as a first-order approximation for the behaviour of a particular queue in the A^-queue 
case, switchover periods now also representing glue and visit periods at other queues. In Section]^ 


the A^-queue case is analyzed. Section presents some conclusions and suggestions for future 
research. 

2 Queue length analysis for the single-queue case 

2.1 Model description 

In this section we consider a single queue Q in isolation. Customers arrive at Q according to a 
Poisson process with rate A. The service times of successive customers are independent, identically 
distributed (i.i.d.) random variables (r.v.), with distribution B{-) and Laplace-Stieltjes transform 
(LST) B{-). A generic service time is denoted by B. After a visit period of the server at Q it 
takes a vacation. Successive vacation lengths are i.i.d. r.v., with S a generic vacation length, with 
distribution S{-) and LST S{-). We make all the usual independence assumptions about interarrival 
times, service times and vacation lengths at the queues. After the server’s vacation, a glue period 
of deterministic (i.e., constant) length begins. Its significance stems from the following assumption. 
Customers who arrive at Q do not receive service immediately. When customers arrive at Q during 
a glue period G, they stick, joining the queue of Q. When they arrive in any other period, they 
immediately leave and retry after a retrial interval which is independent of everything else, and 
which is exponentially distributed with rate v. The glue period is immediately followed by a visit 
period of the server at Q. 

The service discipline at Q is gated: During the visit period at Q, the server serves all "glued" 
customers in that queue, i.e., all customers waiting at the end of the glue period (but none of those 
in orbit, and neither any new arrivals). 

We are interested in the steady-state behaviour of this vacation model with retrials. We hence 
make the assumption that p := AES < 1; it may be verified that this is indeed the condition for 
this steady-state behaviour to exist. 

Some more notation: 

Gn denotes the nth glue period of Q. 

Vn denotes the nth visit period of Q (immediately following the nth glue period). 

Sn denotes the nth vacation of the server (immediately following the nth visit period). 

X„ denotes the number of customers in the system (hence in orbit) at the start of G„. 

Vn denotes the number of customers in the system at the start of Vn- Notice that here we should 
distinguish between those who are queueing and those who are in orbit: We write Vn = Vn'^^ -\-Vn°\ 
where q denotes queueing and o denotes in orbit. 

Finally, 

Zn denotes the number of customers in the system (hence in orbit) at the start of S'„. 

2.2 Queue length analysis at embedded time points 

In this subsection we study the steady-state distributions of the numbers of customers at the 
beginning of (i) glue periods, (ii) visit periods, and (hi) vacation periods. Denote by A a r.v. with 
as distribution the limiting distribution of A„. V and Z are similarly defined, and V = V^'^'> 
the steady-state numbers of customers in queue and in orbit at the beginning of a visit period 
(which coincides with the end of a glue period). In the sequel we shall introduce several generating 
functions, throughout assuming that their parameter \z\ < 1. For conciseness of notation, let 
I3{z) := B{\{1 — z)) and cr(z) := 5(A(1 — z)). Then it is easily seen that 


E[z^] = a{z)E[z^\, 


( 2 . 1 ) 


since X equals Z plus the new arrivals during the vacation; 


since Z equals V'^°^ plus the new arrivals during the services; and 


y(9) y(°) 


( 2 . 2 ) 



(2.3) 


The last equation follows since is the sum of new arrivals during G and retrials who return 
during G; each of the X customers which were in orbit at the beginning of the glue period has a 
probability 1 — e~''^ of returning before the end of that glue period. 


Combining Equations (2.1)-(2.3), and introducing 

fiz) := 




(2.4) 


we obtain the following functional equation for £[ 0 ^]: 

E[z^] = 

Introducing K{z) := and X{z) := E[z^], we have: 


X{z)=K{z)X{f{z)). 


(2.5) 


This is a functional equation that naturally occurs in the study of queueing models which have 
a branching-type structure; see, e.g., |S] and ED- Typically, one may view customers who newly 
arrive into the system during a service as children of the served customer ("branching"), and 
customers who newly arrive into the system during a vacation or glue period as immigrants. Such 
a functional equation may be solved by iteration, giving rise to an infinite product - where the 
jth term in the product typically corresponds to customers who descend from an ancestor of j 
generations before. In this particular case we have after n iterations: 

n 

X(z) = n K{f(^\z))X{f(-+^\z)), (2.6) 

3=0 

where /*'°^(z) := z and f^^\z) := j = 1,2,.... Below we show that this product 

converges for n —)■ cx) iff p < 1, thus proving the following theorem: 

Theorem 1. If p < 1 then the generating funetion X(z) = E[z^] is given by 

00 

X{z) = l[K{f^^\z)). (2.7) 

1=0 


Proof. Equation ( |2.5[ ) is an equation for a branching process with immigration, where the number 
of immigrants has generating function K{z) and the number of children in the branching process has 
generating function f{z). Clearly, K'{1) = AES' -I- \pG < 00 and /'(I) = + (l ~ p < 1, 

if p < 1. The result of the theorem now follows directly from the theory of branching processes 
with immigration (see e.g.. Theorem 1 on page 263 in Athreya and Ney [5]). □ 


i[^ expressions for E,[z^ 


and 


Having obtained an expression for E[z^] in Theorem.^. j_j 

immediately follow from (2.2l and (2.3l. Moments of X may be obtained from Theoremfb but it is 
also straightforward to obtain EA from Equations 


EX = AES -h EZ, (2.8) 

EZ = pEy(«)-hEF(°), (2.9) 

Ey(«) = AG-h (1 - e-'^^)EA, (2.10) 

Ey(o) = e-^‘^EA, (2.11) 


yielding 


Hence 


AES + A^_ 

(l-p)(l-e-‘^G) 


Er(«) = AG-h (1 


AES ApG 

(1 - p)(r- e-G) 


AES -h AG 

1-p 


( 2 . 12 ) 


(2.13) 











(2.14) 


Ey(o) ^ 


EZ = 


ApG + AE5'[p(l - 


(l_p)(l_e-G) • (2-15) 

Notice that the denominators of the above expressions equal 1 — Also notice that it makes 

sense that the denominators contain both the factor 1 — p and the probability 1 — that a 
retrial returns during a glue period. 

In a similar way as the first moments of X, Y^'^\ Y^°'> and Z have been obtained, we can also 
obtain their second moment. For further use we here calculate E[A(A — 1)]: 

K"{1) 


E[A(A- 1)] = 


+ 


(1 - p){l - e-'^G)(i ^(1 _ e-i.G) e-i.G) 

K'{l)[2K'{l){p{l - e-^^) + e-^-^) + (1 - e-'^^)X^EB^] 


(2.16) 


(1 - p)2(l - e-"G)2(i + p(i _ e—G) + e-^'G) 

where A:'(1) = AE5 + XpG and K”{1) = A^ES'^ + 2pX^GES + X^GEB'^ + {XGp)'^. 

Remark 1. Special cases of the above analysis are, e.g.: 

(i) Vacations of length zero. Simply take a{z) = 1 and ES" = 0 in the above formulas. 

(ii) 1 / = oo. Retrials now always return during a glue period. We then have f{z) = /3(z), which 
leads to minor simplifications. 

Remark 2. It seems difficult to handle the case of non-constant glue periods, as it seems to lead 
to a process with complicated dependencies. If G takes a few distinct values Gi,..., G^r with 
different probabilities, then one might still be able to obtain a kind of multinomial generalization 
of the infinite product featuring in TheoremOne would then have several functions fi{z) := 
(1 — e~'''^')f3{z) + e~'^^'z, and all possible combinations of iterations fi{fhifk{- ■ ■ (^)))) arising in 
functions Ki{z) := , i = 1, 2,..., iV. By way of approximation, one might stop 

the iterations after a certain number of terms, the number depending on the speed of convergence 
(hence on 1 — p and on 1 — e”"^®*). 


2.3 Queue length analysis at arbitrary time points 

Having found the generating functions of the number of customers at the beginning of (i) glue 
periods (E[ 2 ;^]), (ii) visit periods (E[z^'’'zJ*°']), and (hi) vacation periods (E[ 2 :'^]), we can also 
obtain the generating function of the number of customers at arbitrary time points. 

Theorem 2. If p < 1, we have the following results: 

a) The joint generating function, Rva{zq, Zo), of the number of customers in the queue and in the 
orbit at an arbitrary time point in a vacation period equals the generating function Rva{zo) of 
the number of customers in orbit at an arbitrary time point in a vacation period and is given by 


Rva{Zq,Zo) = E[z„ ] 


^.l-^(A(l-z„)) 


(2.17) 


X{l-Zo)ES ■ 

b) The joint generating function, Rgi{zq,Zo), of the number of customers in the queue and in the 


orbit at an arbitrary time point in a glue period is given by 

rG 


Rgl i.Zq, Zq ) — 


e-^(i-"’)*E[{(l-e-"*)z5 






(2,18) 


cj The joint generating function, Rm(zq, Zo), of the number of customers in the queue and in the 
orbit at an arbitrary time point in a visit period is given by 


RviiZq, Zo) — 


Zq 

- E[H(A(1 - 

l-B{X{l-Zo)) 


E[r(9)] (z, - R(A(1 - Zo))) 

A(1 - Zo)EB 


(2.19) 











d) The joint generating function, R{zq, Zq), of the number of customers in the queue and in the 
orbit at an arbitrary time point is given by 

R{Zq, Zq') — pRvi (^Zq, Zq') “t“ (1 p') Rgl ^o'} (1 p') Rva{^q ? ^o) • (2.20) 


Proof. 


a) Follows from the fact that during vacation periods there are no customers in the queue, hence 
the right-hand side of (2.17) is independent of Zq, and the fact that the number of customers 
at an arbitrary time point in the orbit is the sum of two independent terms: The number of 
customers at the beginning of the vacation period and the number that arrived during the past 
part of the vacation period. The generating function of the latter is given by 


1 - 5(A(1 - zo)) 
A(1 - Zo)ES 


b) Follows from the fact that if the past part of the glue period is equal to t, the generating 
function of the number of new arrivals in the queue during this period is equal to 

and each customer present in the orbit at the beginning of the glue period is, independent of 
the others, still in orbit with probability e~'^^ and has moved to the queue with probability 
1 -e"^*. 

c) During an arbitrary point in time in a visit period the number of customers in the system 
consists of two parts: 

• the number of customers in the system at the beginning of the service time of the customer 
currently in service, leading to the term (see Remark]^ below): 


z, - e[b(a(i - zo))^‘^r^°’]) 

E[y(9)] (^Zq - B(A(1 - Zo))) 


( 2 . 21 ) 


• the number of customers that arrived during the past part of the service of the customer 
currently in service, leading to the term 


1 - .B(A(1 - Zq)) 

A(1 — Zo)Ei3 

d) Follows from the fact that the fraction of time the server is visiting Q is equal to p, and if 
the server is not visiting Q, with probability ES'/(G + ES”) the server is on vacation and with 
probability G/(G -I- ES") the system is in a glue phase. 

□ 


Remark 3. A straightforward way to prove (2.21 1 is to condition on the number of customers, say, 
j, in queue at the end of a glue period, and to average the number of customers in queue and 
in orbit over the beginnings of the j services. A more elegant proof of Formula (2.21) uses the 
observation that typically in vacation and polling systems each time a visit beginning or a service 
completion occurs, this coincides with a service beginning or a visit completion. This observation 
yields (see, e.g., Boxma, Kella and Kosinski [^) 


7V6(z,, Zo) -k Sc{Zq, Zo) = Sb{Zq, Zo) + jVc{Zq, Zo). 

Here, Vb{zq, Zq) and Vc(Zg, Zq) are the joint generating functions of the number of customers in the 
queue and in the orbit at visit beginnings and visit completions, respectively. Similarly, Sb{zq,Zo) 
and Sc{zq, Zo) are the joint generating functions of the number of customers in the queue and in 
the orbit at service beginnings and service completions, respectively. Finally, 7 is the reciprocal of 
the mean number of customers served per visit. Clearly, 

7 = > V6(z„Zo) = E[z^'‘'’zj'°’], Vc(z„Zo) =E[H(A(l-Zo))^'''’zj'°’], 











and 


^o) — 


^b\^qj^o) 


B{X{l-Zo)), 


which yields that Sb{zq,Zo) is given by (2.211. 


From Theorem we now can obtain the steady-state mean number of customers in the system 
at arbitrary time points in vacation periods in glue periods (E[i?g;]), in visit periods 

(E[l?„i]) and in arbitrary periods (E[i?]). These are given by 


E[R,4 = E[Z] + A|§| 
E[i?g,] =E[X] + A« 


2 ’ 


E[i?„,] = l + A|gj + 

E[i?] = pE[ii„i] -I- (1 - p)(3+E5E[i?g/] -I- (1 - p) 


Efrc®)] 2E[y(9)] ■ 

G TirrD 1 I /-I ES 


( 2 . 22 ) 


Remark that the quantities and — 1)] can be obtained using (2.31: 

E[y(9)y(o)] = AGe-^‘^E[X] + (l - e-'^‘^E[X(A: - 1)], 

E[r(9)(y(9) _ 1)] = (ag) 2 -f (1 -e-‘'‘^)^E[X(X- 1)] -f 2AG(l -e-'^‘^)E[X]. 


By combining these relations with ( |2.22[ ), ( |2.8[ |, ( |2.12 1 and ( 2.16[ ), we obtain - after tedious 
calculations - the following relatively simple expression for the mean number of customers E[i?]: 




A2E[B=| AE[(G + S)=1 , ApE|G + S] , 
'^1"' * 20^ + 2EIG + S] + 1-f, + 


(2.23) 


which we rewrite for later purposes as 

E|fi] =, + » + A JH|L ss + A JHILg + A 


2(1-p) G-f E)^] 2E[S'] 
Ap(G + E[5]) 


'G-f E)^] 


G G 
G + E)^] '2 


1- p 


A(pG + E[5]) 


(l-p)(l-e-G)- 


(2.24) 


Remark 4- (i)- It should be noticed that the first two terms in the right-hand side of (2.231 together 


represent the mean number of customers in the ordinary MjGjl queue, without vacations and 
glue periods. The third term represents the mean number of arrivals during the residual part of 
a vacation plus glue period. The fourth term can be interpreted as the mean number of arrivals 
during a visit period of the server (since the mean length of one cycle, i.e., visit plus vacation plus 
glue period, is via a balance argument seen to equal E[G] = while a mean visit period 

equals pE[G]). The fifth term is the only term involving the retrial rate v. In particular, that term 
disappears when ;/ — )• cxd. In the latter case, our model reduces to an MjGjl queue with gated 
vacations, with vacation lengths G -I- S'. The resulting expression for the mean number of customers 
coincides with formula (5.23) of |53] (see also formula (3.2.6) of |27|1. 

(ii). A more interesting limiting operation is to simultaneously let z/ —> oo and G ), 0, such that vG 
remains constant. The resulting model is an MjGjl queue with vacations and binomially gated 
service; see, e.g.. Levy |15| . Here, each customer who is present at the end of a vacation, will be 
served in the next visit period with probability p = 1 — . In this case, the mean number of 

customers in the system is given by 


A“E[B2] AE[S2] ApE[Sl , AE|S](l-p) 

+ 2(wj + ^ + ^iT-yr- 


(2.25) 


This formula coincides with the results obtained in [15] (see e.g., formula (7.1) with A7 = 1 in [15] 
for the mean sojourn time of customers in this model). Observe that our formula, after application 





































of Little’s formula, does not fully agree with the mean delay expression (5.50b) in |53| and with a 
similar formula on p. 90 of m- Those mean delay expressions for the binomial gated model seem 
to refer to the case where customers who are not served during a visit (w.p. 1 — p) are lost; hence 
factors like — in those mean delay expressions. 


iii). Formula (2.241 immediately shows how the mean number of customers behaves for very small 


and for very large values of the glue period length G: 


and 


E[R] 


E[i?] 


AE[S'] 

GH^-P)’ 


A(1 + p) 

2(1-p) 


G, 


GiO, 

(2.26) 

G —>■ oo. 

(2.27) 


We observe that the equations (2.261 and (2.271 do not involve E[i3^] and E)^^]. Hence the mean 
number of customers in the system is almost invariant to the variance of the service times and 
vacation times when G is either very small or very large. In Subsection |2.4| we explore the effect of 
G on E[i?] more deeply. 


2.4 Optimizing the length of the glue period 

One of the main reasons for studying mathematical models of optical networks is to improve the 
performance of the system. In this model the length of the glue period is an important system 
design parameter. The results of the previous subsections can, e.g., be used to determine the value 
of G which minimizes the mean number of customers in the system at any arbitrary time point. 
The mean sojourn time of an arbitrary customer follows from Little’s formula. Therefore we can 
find the value of G which minimizes the mean sojourn time of an arbitrary customer. 

Let us first present a sample of numerical results that we obtained for E[i?] as a function of 
G. We consider four cases: (i) the service time distribution and vacation time distribution are 
exponential, (ii) the service time distribution is exponential and the vacation time is constant, (iii) 
the service time is constant and the vacation time distribution is exponential and (iv) the service 
time and vacation time are constant. In Figwe take A = 0.9, z/ = 0.5, and ES' = 10, we plot G vs 
E[i?], for EB = 1 and 1.1. Note that EB = 1.1 corresponds to the heavy traffic case p = 0.99. 
From the examples in Fig we observe the following results: 

— The mean number of customers at any arbitrary point seems to be convex w.r.t. glue period 
length, i.e., there exists a glue period Gmin where the system has minimum mean number of 
customers E[Rrain\ and hence minimum mean sojourn time. 

— For a large G, E[i?] increases linearly with G (as confirmed 

— For a very small G, E[i?] behaves like 1/G (as confirmed by 

— For different service time and vacation time distributions the mean number of customer is 
changed but the difference in Gmin is insignificantly small. 

— As p —>■ 1 , i.e. the system is in heavy traffic regime, the value of E[i?] becomes insensitive to 
the distribution of S. 


b y d^ ). 

(pfT 


Indeed, if G is very small, the number of customers joining the queue in each glue period is 
very small and thereby the efficiency of the system is decreased. On the other hand, a large G 
means the system stays in the glue period for a long time and this decreases the efficiency of the 
system. Hence it is logical to have a Gmin which optimizes the system. 

We will now prove that E[R] is indeed convex in G. Twice differentiating the expression for 
E[i?] in Equation ( 2.23[ ) w.r.t. G gives 


dG"^ 


E[i7] 


A 

E[G + S] 


E[(G + 5)2] 
E[G + S']2 


^ (l-p)(l-e—G)2 


i^{pG + E[S]) 


(1 — 


-2p . 
(2.28) 



























(a) EB = 1 



Fig. 1: Mean number of customers vs length of glue period 


Clearly, the first term in the right-hand side of (2.28) is nonnegative. Let 


Q(G) :=KpG + E[S]) 


(1 




(1 — 


-2p. 


We can see Q{G) —oo as G —0 or G —> oo. Let Q{G) attain its minimum at G = g. Hence 
calculating the first derivative of Q{G) and equating it to zero, at G = g, gives 


P5 + E[^] = 


P 


2ve-'^3 


(l-e-2"s) 


Therefore 


Q(5) = |[e"®-2 + e-^®]>0. 


We observe that the minimum value of Q{G) is always nonnegative. Since both terms of (2.28) are 
nonnegative, -^^E[R] > 0. 

Hence E[i?], the mean number of customers at an arbitrary point of time in the system, is 
convex in G. So the system can improve the quality of service by setting an optimal value of G for 
the fixed glue period. 

In Tablewe analyze the behaviour of Gmin and E[Rmin\ as we increase £[5”] for an exponential 
distribution B{-) with EH = 1, arrival rate A = 0.5 and retrial rate v = 0.5. 

Table [2 suggests that, in the case under consideration. 


— Gmin and E[Rmin\ increase when £[5] increases. 

“ Gmin and E[Rmin\ increase when the variance of S becomes larger. 

— When EH approaches 0, Gmin also approaches 0. When there is no customer in the queue, the 
system will then have a series of very short glue periods, and when a customer arrives or returns 
from orbit, it can almost instantaneously be taken into se rvice. In this case, the system reduces 
to an ordinary M/G/1 retrial queue; indeed. Formula (2.23) reduces to E[i?] = p + 


which is in agreement with Formula (1.37) of |11| . 

























Table 1: Optimal glue period length and corresponding minimum mean number of customers 


ES 

Grain for 
exponential S{-) 

Grain foj' 
constant S{-) 

E[R„,in] for 
exponential S{-) 

E[Amjnl for 
constant S{-) 

0 

t, t —0 

o 

T 

~ 2 

~ 2 

0.1 

0.608 

0.606 

2.476 

2.473 

0.5 

1.334 

1.320 

3.419 

3.385 

1 

1.846 

1.801 

4.259 

4.170 

5 

3.694 

3.508 

9.276 

8.549 

10 

4.785 

4.495 

14.778 

13.070 

50 

7.783 

7.225 

56.027 

45.156 

100 

9.167 

8.521 

106.606 

83.634 


3 Queue length analysis for the N-queue case 

3.1 Model description 


In this section we consider a one-server polling model with multiple queues, Qi, i = 1, • • • ,N. 
Customers arrive at Qi according to a Poisson process with rate A^; they are called type-z customers, 
z = 1, • • • , N. The service times at Qi are i.i.d. r.v., with Bi denoting a generic service time, with 
distribution Bi{-) and LST 13i{-), i = 1, - ■ ■ ,N. The server follows cyclic polling, thus after a visit 
of Qi, it switches to Qi+i- Successive switchover times from Qi to Qi+i are i.i.d. r.v., with Si a 
generic switchover time, with distribution Si{-) and LST S'i(-), z = I,-- - ,N. We make all the 
usual independence assumptions about interarrival times, service times and switchover times at 
the queues. After a switch of the server to Qi, there first is a deterministic (i.e., constant) glue 
period Gi, before the visit of the server at Qi begins, z = 1, • • • , TV. As in the one-queue case, the 
significance of the glue period stems from the following assumption. Customers who arrive at Qi 
do not receive service immediately. When customers arrive at Qi during a glue period Gi, they 
stick, joining the queue of Qi. When they arrive in any other period, they immediately leave and 
retry after a retrial interval which is independent of everything else, and which is exponentially 
distributed with rate Vi, i = 1, ■ • • , TV. 

The service discipline at all queues is gated: During the visit period at Qi, the server serves all 
"glued" customers in that queue, i.e., all type-z customers waiting at the end of the glue period - 
but none of those in orbit, and neither any new arrivals. 

We are interested in the steady-state behaviour of this polling model with retrials. We hence 
assume that the stability condition Pi < 1 bolds, where pi := XiEBi. 

Some more notation: 

Gni denotes the nth glue period of Qi. 

Vni denotes the rzth visit period of Qi. 

Sni denotes the nth switch period out ol Qi, i = 1, • • • , TV. 

, X ^2 i ■ ■ ■ : ^nlf) denotes the vector of numbers of customers of type 1 to type TV in the system 
(hence in orbit) at the start of Gni, z = 1, • • • , TV. 

^^ 2 ^ ■ ■ ■ ) denotes the vector of numbers of customers of type 1 to type TV in the system 
at the start of Vni, z = 1, • • • , TV. We distinguish between those who are queueing in Qi a-nd those 
who are in orbit for Qp. We write -I- YnT\ where q denotes queueing 

and o denotes in orbit. 

Finally, 

z\^ 2 , • • • , Z^]q) denotes the vector of numbers of customers of type 1 to type TV in the system 
(hence in orbit) at the start of Sni, z = 1, • • • , TV. 




3.2 Queue length analysis 


In this subsection we study the steady-state joint distribution of the numbers of customers in 
the system at beginnings of glue periods. This will also immediately yield the steady-state joint 
distributions of the numbers of customers in the system at the beginnings of visit periods and 
of switch periods. We follow a similar generating function approach as in the one-queue case, 
throughout making the following assumption regarding the parameters of the generating functions: 
\zi\ < 1, \ziq\ < 1, \zio\ < 1. Observe that the generating function of the vector of numbers of 
arrivals at Qi to Qjv during a type-i service time Bi is Z 2 , • • • , Zn) := Aj(l — zj)). 

Similarly, the generating function of the vector of numbers of arrivals at Qi to Qn during a type-j 
switchover time Si is ai{zi, Z 2 , - ‘ ‘ ) zn) := Si{J2f=i ~ Since the server serves Qi+i after 




Qi we can successively express (in terms of generating functions) ■ ■ ■ 

I yit) y{i) yii) \ ( y(i) y(i) x ■ . vhl) \r{io} t QTirl 

! -^n2 1 ■ ■ ■ J i ^n2 J ’ ’ ’ j Onl ; ^ n2 t ^ ni ^ ^ ni J ’ ’ ’ i ^nNJ’ 
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,y«)into {x^lx^l- 


X^]q); etc. Denote by 

the vector with as distribution the limiting distribution of ^^ 2 ^ ’ ’ ’ ? ’ ’ ’ > 


I ) 


and similarly introduce ■ 

for z = 1, • • • ,N. We have: 


and with Y^^^> = + Y^ 
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Furthermore, 
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. (3.3) 


It follows from (3.1 1 , (3.2) and (3.3), with 


fi{zi,Z2,-■ ■ ,zn) ■■= {I-0 '''^Qj3i[zi,z2,-" ,Z]si) + e ‘''^'Zi, 
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(3.4) 


Since the server moves to Qi after Qn, substituting i = in (3.41, we have 
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From (3.4) we have 


E 


n 

Kj9^N 


X 


(N) 


[hN{z)] 


X 


(N) 


' N-2 


= a<-^-^\z) I Yl 

J=-L 


X E 


'N-2 

n 

,i=i 


AN-l) 


[hN-l{z)f--^ [h^{z)f- 
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From (3.5) and (3.6) we have 
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where 


2-1 N 


i=i 

By recursively substituting as above we get 

N N 
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r(i) 




j=i+l 




^(1) 


r(i) 


i=l 


i=l 


Equation (3.7) can be divided into three factors, representing the switchover period, glue period 
and visit period respectively. The first factor, for a particular value of i, represents the arrivals 
during the switchover time after the visit of Qi. The second factor represents the arrivals during the 
glue period before a visit of Qi. It is further divided into three generating functions. First are the 
arrivals of type j < i; these don’t have any further effect on the system. Then the arrivals of type 
i, these are served during the following visit and produce new children (i.e., arrivals during their 
service) of each type. Finally those of type j > i which may or may not be served in future visits 
and if served produce new children of each type. These two factors are taken for all i = 1, • • • , N. 
The third factor represents the descendants (arrivals during services, arrivals during services of 
customers who arrived during services, etc.) of (-^^^\ • 

Consider 


)■ 


N 

n 




( 1 ) 


X{z)=E 

with an obvious definition of K{z), we can rewrite ( |3.7[ ) into 

X{z) = K{z)X{h{z)) 


(3.8) 


where 


Define, for alH = 1, • • • ,N, 


h{z) := (hi(z),-- - ,hN{z)). 


h'f’iz)=Zi, h'f\z) = '^\z),h ^2 ^\z),---,h^N 

Theorem 3. If^iPi < 1, then the generating function X{z) is given by 

OO 

X{z) = K{h^f^\z),h^^\z),--- ,h^^\z)). 

ra—Q 


(3.9) 


Proof. Equation (3.9) follows from (3.8) by iteration. We still need to prove that the infinite 
product converges if 'YhiPi < 1- Equation (3.8) is an equation for a multi-type branching process 
with immigration, where the number of immigrants of different types has generating function K{z) 
and the number of children of different types of a type i individual in the branching process has 
generating function hi{z), i = 1, • • • ,N. An important role in the analysis of such a process is 
played by the mean matrix M of the branching process, 


mil • • -miN 


M = 


rriNi ■ ■ -rriNf^ 


(3.10) 


where m^ represents the mean number of children of type j of a type i individual. The elements of 
the matrix M are the same as given in Section 5 of Resing m, which is 


N 

m^j = fij • l[j < i] + ^ fiknikj, 
k=i+l 


(3.11) 

















where = |^( 1 , 1 , • • • , 1 ) and /„ = |^( 1 ,1, • • • , 1 ). 

We observe that the equation for rriij is the sum of two terms. First the children of type j < i, 
who do not affect the system in the future. Next the children of type j produced by the children of 
type fc > i in the subsequent visits. 


The theory of multi-type branching processes with immigration (see Quine |20| and Resing El) 
now states that if (i) the expected total number of immigrants in a generation is finite and (ii) the 
maximal eigenvalue Xmax of the mean matrix M satisfies Xmax < Ij then the generating function 
of the steady state distribution of the process is given by (3.9). To complete the proof of Theorem 
1 ^ we shall now verify (i) and (ii). 


Ad (i): The expected total number of immigrants in a generation is 
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E5, 
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m. 


■ij 


E(G,+E5,)e 

i=i 


^E5, + ^ G, . 

j=i+l J 

(3.12) 




Since the above equation is a finite sum/product of finite terms it is indeed finite. 

Here, the term Ai(Gi(^^-mij)) corresponds to the type 1 customers arriving during the 
glue period of Qi and their subsequent children of all types. The term 

corresponds to the type 1 customers arriving during the glue periods of Qj, j = 2, • • • ,7V, and 
switchover periods after Qj, j = - ■ ,N. These customers arrive after the visit of Qi and hence 

do not get served or produce children. The term A 2 (((Gi -b ES'i)(l — -b G 2 )(X)j ^ 2 j) + 

(Gi -bEiS'i)e“'^^‘^^) corresponds to the type 2 customers arriving during the glue period of Qi, Q 2 , 
the switchover period after Qi and their subsequent children. The term A 2 (X]j =2 Si =3 ) 

corresponds to the type 2 customers arriving during the glue periods of Qj, j = 3, • • • ,7V, and 
switchover periods after Qj , j = 2, • • • ,7V. These customers do not produce any children. Similarly 

the term AAr(((EfE(Gj + E^,))(l - + Gn^Ej^Nj) + EjJiiGj + 

corresponds to the type TV customers arriving during the glue period of Qi, • • • , Qat, the switchover 
periods after Qi, - ■ ■ , Qn-i their subsequent children. The term AjvES'tv corresponds to the 
type TV customers arriving during the switchover period after Qn, which do not produce any 
children. 

Ad (ii): Define the matrix 


/e-i'iGi + (1 - p-^ 

(1 _e-'^2G2) A 1 EH 2 

(l_e-'2iGi) A 2 ER 1 
g-i22G2 • 

(1 - AatERi \ 

Xn^B2 

(1 - AiERat 

(1 _ e-'^wGjv) A2ERAr 

■ + [l- pn) 


(3.13) 


where the elements hij of the matrix H represent the mean number of type j customers that 
replace a type i customer during a visit period of Qi (either new arrivals if the customer is served, 





or the customer itself if it is not served). We have that 
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/ebA 

EB2 


\Ei?Ary 



Q-U 2 G 2 



{Y,P^) 

EBi \ 


/ERi\ 


EB2 


EB2 
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< 


Pj) 

EBn j 


\EBm) 


(3.14) 


if and only if pj < 1. Using this result and following the same line of proof as in Section 5 of 
Resing |21) . we can show that the stability condition Yj Pj < 1 implies that also the maximal 
eigenvalue Xmax of the mean matrix M satisfies Xmax < 1- This concludes the proof. □ 


We can now obtain the moments, either from (3.91 or in a similar way as in Section 


2.2 in terms of EXj* 

and Exf ^ : 



EX 

=X,ES,+EZf\ 

When j ^ i, 

EX]*^ 

= XjEBiEvY^ + Ey/*\ 

else 

Exf^ 

= a,er,ey/*«^ + ey/*°^ 

Further 

ey/*^ = 

= AjG,+Exf\ 


eyY^ = 

Eybo) ^ 

= A,G, + (l-e-"-^-)Ex| 
e"'^*®‘EXf\ 


From the above equations we get, when j ^ i ■ 


= XjES, + Aj(l + p,)G, + XjEB,{l - e-"-^’)Ewf ^ + Exf, 


and 


Exf+^) = X,ES, + X,p,G, + (p,(l - e-"’^-) + e-"'^-)Exf. 

Using flow balance arguments (mean number of customers of type i served per cycle equals 
mean number of type i customers arriving per cycle) and the obvious fact that the mean cycle 
time equals EC := + Gi)/{1 — p) , we obtain 




(3.15) 


We can also use a similar argument for mean number of type i customers leaving the orbit, 
(1 — e“'^’^’)EX|*\ to equal the mean number of type i customers entering it, Ai(EC — Gi), per 
cycle, yielding 


EX.«> = A. ^ [Ej(ES, +G,) _ 


1 .G L 1 (3-16) 

1 — 1 — P 

A - V A -k Hq) (i) 

We can observe that (3.151 and (3.16) satisfy the above relation between EUl and EX) . Further 

'-' '-' (i) (i) (i) 

for each i, cyclically substituting we get all EXj and therefore EY^ and EX) . 

The second moments of xj^^ and the various terms E[x[^X^*^] can be obtained by solving a 
set of equations which is derived by twice differentiating (3.8) w.r.t. Zj and Zk, j,k = 1 , - ■ ■ ,N, 




















and calculating the value at z 


(iV 

(I,!,-- - ,1). Since the system is cyclic, once we obtain 


(i)2 

j = 1, • • • ,N,we can similarly obtain EXj , z = 1, • • • , iV, by changing indices. It is not difficult to 
develop an efficient procedure for determining higher moments in polling systems with a branching 
discipline, cf. m- 


3.3 Queue length analysis at arbitrary time points 


In the previous section we have given the procedure for finding the distribution of the number of cus- 


r{i) 


tomers at the beginning of (i) glue periods (E[J|^ Zj ^ ]), (ii) visit periods (E[ ( Ylj^i ^ ) ^iq 


Y, 


(i) 


Ai) 


and (iii) switchover periods (E[J^^- Zj ^ ]), for z = 1, • • • , N. Similar to the vacation model, we now 
obtain the generating functions of the numbers of customers in queue and in the orbit, at all the 
stations, at arbitrary time points. 


Theorem 4. If Pz < 1 [zq, Zo) '■= (zi^, Zio, • • • , zjvq, zato), we have the following results: 

a) The joint generating function, R^fw{zq, Zg), of the numbers of customers in the queue and in 
the orbit at an arbitrary time point in a switchover period after Qi equals the joint generating 
function, Ri^w{zo), of the numbers of customers in orbit at an arbitrary time point in a switchover 
period after Qi and is given by 


'l-FfTT ^jo)) 

"E^ 




(3.17) 


b) The joint generating function, R^^^ {zq, Zo), of the numbers of customers in the queue and in 
the orbit at an arbitrary time point in a glue period of Qi is given by 


Rfii.zq,zo) = 










(3.18) 

c) The joint generating function, R^^}[zq, Zo), of the numbers of customers in the queue and in 
the orbit at an arbitrary time point in a visit period of Qi is given by 


Rfl{Zq,Zo) = 


Ziq I Hzlj 
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(3.19) 


d) The joint generating function, R{zq, Zo), of the numbers of customers in the queue and in the 
orbit at an arbitrary time point is given by 

R{Zq, Zo) = (^PiR\^}{Zq, Zo) + (1 — Pi) ^^q, Zo) + (1 — Pi) (Gj+ES„-) 42 J ^o)^ • 

(3.20) 


Proof. The proof follows the same lines as the proof of Theorem in particular for parts a and d. 
We restrict ourselves here to an outline of the proof of parts b and c. 







b) Follows from the fact that if the past part of the glue period is equal to t, the generating 

function of the number of new arrivals of type i in the queue during this period is equal to 
g-Ai(i-zi,)t each type i customer present in the orbit at the beginning of the glue period 
is, independent of the others, still in orbit with probability and has moved to the queue 
with probability 1 — Further the generating function of the number of new arrivals of 

any type j ^ i in the queue during this period is equal to . 

c) During an arbitrary point in time in a visit period the number of customers in the system 
consists of two parts: 

• the number of customers in the system at the beginning of the service time of the customer 
currently in service, leading to the term 
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j^t ^jo 
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(see Remark]^. 

• the number of customers that arrived during the past part of the service of the customer 
currently in service, leading to the term 


1 ho)) 

(Ej ho)) E-Sz 


From Theorem we now can obtain the steady-state mean number of customers in the system 
at arbitrary time points in switchover periods (E[i?i2]) after Qi, in glue periods (E[i?^j^]) and in 

visit periods (E[_r|,*^]) of Qi , for i = 1, • • • , TV, and at any arbitrary time point (E[i?]). These are 
given by 


3 

3 

E[^ii^] = 1 + (E ^ 

j ' jh 

E[i?] = E {p^Rvh + (1 - P)OT^E[i?«] + (1 - • 




(3.21) 


The mean number of type k customers in the system at arbitrary time points in a switchover period 
after Qi and a glue period before Qi are given by the values of the fc-th term of the sums in the 
formulas of E[i?i*2] and E[i?g)^]. The mean number of type i customers in the system at arbitrary 
time points in a visit period of Qi is given by 


*2E[Bi] Efyl'’'] 2 E[y1‘'’^] 


The quantities and 


1)] can be obtained using (3.31. 


Remark 5. Using a similar approach as presented in m for a polling system with binomial-gated 
service, we can also obtain the following expression for E[iii], the steady-state mean number of 
type-i customers in the system at arbitrary time points. 


E[i?i] — Pi + 
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Efy/'"’] 
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(3.22) 

















which, after summing over z, leads to the alternative formula 
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(3.23) 


Remark 6. From (3.23), we can derive an explicit expression for the mean number of customers 


in the system in the case of a completely symmetric system (A^ = A/iV, pi = p/N, EBi = EB, 
EBf = EB^, ES'i = £5, ESf = ES'^, Gi = G,v, = v). In this case we get 


E[i?] — p + 


X^E[B'^] XN{G + ES) , AVar(S') 

2(1 -P)^ 2 ^ 


{N + l)Xp{G + ES) Xe-^^ 

2(1 - p) 1 - 


2{G + ES) 
N{G + ES') 


1 - p 


-G 


(3.24) 


Remark 7. In [7] the following so-called pseudo conservation law - an explicit expression for 
E PjE[hFi]) with E[Wi] the mean waiting time of a customer of type i until the start of its service - 
has been proven for a large class of polling systems, which also contains the present model: 


^ P.E[1F,] = p%^^M + 


2(1-p) 


2 EE,(^, + G 0 ] 


2(1-p) 




(3.25) 

where Eii^i + Gi) is the sum of all the idle periods of the server and Fi is the work left in Qi at 
the start of a switchover from Qi. Hence, E[F!i] = E[z|*^]E[i?i]. Other than E[Fi], the expression is 
independent of the service discipline. Using a fairly straightforward mean value analysis we obtain 


E[F,] = p^E[C] 
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(£[ 0 ] - GO. 


(3.26) 


From Equations (3.25) and (3.26), we obtain the following pseudo conservation law: 
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(3.27) 


The use of this pseudo conservation law seems to be the easiest way to derive (3.24). Another 
useful aspect of this pseudo conservation law is that it allows us to study the effect of the length of 
the glue period. 

In a more general (not necessarily optical switch related) setting, referred to in the Introduction, 
the glue period may represent the only opportunity to make a reservation for service. The glue or 
reservation period now is the last part of the switchover period to Qp, one could view S* := -I- Gi 


as the total switchover period into Qi and Gi as the reservation period. Formula (3.27) allows us 


to study how the mean waiting times or mean queue lengths are affected by having only a brief 
reservation period, instead of being able to make a reservation at any time (which is the classical 
gated polling system). The only difference in ^ piE[Wi] between our reservation model and the 
classical gated polling model is the last term. If ViGi is very small, that last term will be dominant. 
If it is not very small whereas, e.g., the second moments of the service times are large, then the 
extra term is relatively small - and the advantage for the system operator of having to offer only 
very limited reservation opportunities may outweigh the fact that waiting times and queue lengths 
become slightly larger. 


































3.4 Numerical example 


In this subsection we present some numerical results for the 2-queue case. We have numerically 
evaluated the expressions for E[ii 2 ] and E[i?] using Eq. (3.211. We have also verified that 

Eq. (3.221 gives exactly the same E[i?i] values, and that the pseudo conservation law (3.27) 
is satisfied in the numerical examples. In the numerical example in this section we consider a 
model with two stations. We have chosen Ai = 2, A 2 = 1, = 1 and vi = 0.2. The switchover 

times are deterministic with E[S'i] = 0.5 and E[S' 2 ] = 1. The service times are exponential with 

E[Bi] = E[B2] = 0 . 2 . 

First of all we look at the mean number of type-1 customers in the system, E[i?i], the mean 
number of type-2 customers in the system, E[i? 2 ], and the mean total number of customers in 


the system, E[i?] = E[i?i] 
Gi = G 2 = G (see Fig[^. 


E[i? 2 ], if we vary the length of the glue period G in the case that 



Fig. 2: Mean number of customers vs length of glue period in case Gi = G 2 = G. 


Next we look in Fig at the case where the glue periods Gi and G 2 are unequal. In particular 
we look at the mean number of customers in the system of type-1, type-2 and in total if we vary 
one of the glue periods while keeping the other glue period constant. 

The plots suggest that there is again a unique value of the length of the glue period that 
minimizes the mean total number of customers in the system. They also show that the length of 
the glue period that is optimal for E[i?] is different from the one that is optimal for E[i?i] and the 
one that is optimal for E[i? 2 ]- 


4 Conclusions and suggestions for future research 

In this paper we have studied vacation queues and iV-queue polling models with the gated service 
discipline and with retrials. Motivated by optical communications, we have introduced a glue period 
just before a server visit; during such a glue period, new customers and retrials "stick" instead 
of immediately going into orbit. For both the vacation queue and the N-queue polling model, we 
have derived steady-state queue length distributions at an arbitrary epoch and at various specific 
epochs. This was accomplished by establishing a relation to branching processes. 

Jointly with B. Kim and J. Kim, in [T] we are considering the model with non-constant glue 
periods. We derive E[i?] using MVA (mean value analysis) for the vacation model. It gives us 
the result for general glue period distributions without calculating the generating functions at 
different epochs. We show that a workload decomposition and pseudo conservation law, as discussed 


































(a) varying Gi, with G 2 = 2 


(b) varying G 2 , with Gi = 2 


Fig. 3: Mean number of customers vs length of glue period in case Gi G 2 . 


in [7] and Remark can be derived for these variants and generalizations, and they may be 
exploited for analysis and optimization purposes. We shall then also try to explore the following 
observation: One may view our iV-queue model as a polling model with a new variant of binomial 
gated, with adaptive probability pi of serving a customer at a visit of Qb Pi = 1 when the customer 
arrived in the preceding glue period, and = 1 — otherwise. We would also like to explore 

the possibility to study the heavy traffic behavior of these models via the relation to branching 
processes, cf. [15] . 

From a more applied perspective we are looking into systems with multiple customer classes 
where one class has priority over another. This would help to incorporate the real life scenario 
where some type of data packets, like video buffering, should have as low latency as possible, 
whereas others, like a file transfer, can be delayed a bit longer. 

Finally, we would like to point out an important advantage of optical fibre: the wavelength of 
light. A fibre-based network node may thus route incoming packets not only by switching in the 
time-domain, but also by wavelength division multiplexing. In queueing terms, this gives rise to 
multiserver polling models, each server representing a wavelength. We refer to [2] for the stability 
analysis of multiserver polling models, and to |5] for a mean field approximation of large passive 
optical networks. Therefore we would like to study multiserver polling models with the additional 
features of retrials and glue periods. 
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